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Abstract 

We consider the non-radial motions, originating in the outskirts of 
clusters of galaxies and we show how it may reduce the discrepancy 
between Cold Dark Matter (CDM) predicted X-ray temperature dis- 
tribution function of clusters of galaxies and that really observed, so 
resolving the problem of the X-ray clusters abundance over-production 
predicted by the CDM model. We construct the X-ray temperature dis- 
tribution function using Press-Schechter's (1974) theory and Evrard's 
(1990) prescriptions for the mass-temperature relation, taking also ac- 
count of the non-radial motions originating from the gravitational in- 
teraction of the quadrupole moment of the protocluster. Then, we 
compare our X-ray temperature distribution function with the tidal 
field of the matter of the neighboring protostructures. We find that 
the model produces a reasonable cluster temperature distribution. 

keywords: cosmology: theory-large scale structure of Universe - galax- 
ies: formation 



1 Introduction 



In the last years the CDM model has contributed to obtain a better un- 
derstanding of the origin and evolution of the large scale structure in the 
Universe (White et al. 1987; Frenk et al. 1988; Efstathiou 1990; Ostriker 
1993). The principal assumptions of the SCDM (standard CDM) model are: 
a) a flat Universe dominated by weakly interacting elementary particles hav- 
ing low velocity dispersion at early times; b) Barionic content determined 
by the standard big bang nucleosynthesis model (Kernan & Sarkar 1996; 
Steigman 1996; Olive 1997; Dolgov 1997); c) critical matter density; d) ex- 
pansion rate given by h = 0.5; e) a scale invariant and adiabatic spectrum 
with a spectral index, n = 1; f) the condition required by observations, that 
the fluctuations in galaxy distribution, (8p/p) g , are larger than the fluctua- 
tions in the mass distribution, {5p/ p) p by a factor b > 1. As pointed out by 
Davis et al. (1985), without this last assumption, a unbiased CDM (that is 
6 = 0) presents of several problems: pairwise velocity dispersion larger than 
the observed one, galaxy correlation function steeper than observed (Liddle 
& Lyth 1993 and Strauss Sz Willick 1995). The remedy to these problems 
is the concept of biasing (Kaiser 1984), i.e. that galaxies are more strongly 
clustered than the mass distribution from which they originated. The phys- 
ical origin of such biasing is not yet clear even if several mechanisms have 
been proposed (Rees 1985; Dekel & Silk 1986; Dekel & Rees 1987; Carlberg 
1991; Cen & Ostriker 1992; Bower et al. 1993; Silk k Wyse 1993). 
Although at his appearance the standard form of CDM was very success- 
ful in describing the observed structures in the Universe (galaxy cluster- 
ing statistics, structure formation epochs, peculiar velocity flows) (Peebles, 
1982; Blumenthal et al. 1984; Bardeen et al. 1986; White et al. 1987; Frenk 
et al. 1988; Efstathiou 1990) recent measurements have shown several de- 
ficiencies of the model, at least if any bias of the distribution of galaxies 
relative to the mass is constant with scale. Some of the most difficult prob- 
lems to reconcile with the theory are the magnitude of the dipole of the 
angular distribution of optically selected galaxies (Lahav et al. 1988; Kaiser 
& Lahav 1989), the possible observations of clusters of galaxies with high 
velocity dispersion at z > 0.5 (Evrard 1989), the strong clustering of rich 
clusters of galaxies, £ C c( r ) — (r /25h~ 1 Mpc)~ 2 , far in excess of CDM predic- 
tions (Bahcall & Soneira 1983), the X-ray temperature distribution function 
of clusters, over-producing the observed cluster abundances (Bartlett & Silk 
1993), the conflict between the normalisation of the spectrum of the pertur- 
bation which is required by different types of observations. Normalization 
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obtained from COBE data (Smoot et al. 1992) on scales of the order of 
10 3 Mpc requires ag = 0.95 ± 0.2. Normalisation on scales 10 to 50Mpc ob- 
tained from QDOT (Kaiser et al. 1991) and POTENT (Dekel et al. 1992) 
requires that ag is in the range 0.7 -j- 1.1, which is compatible with COBE 
normalisation while the observations of the pairwise velocity dispersion of 
galaxies on scales r < 3Mpc seem to require ag < 0.5. Another problem of 
CDM model is the incorrect scale dependence of the galaxy correlation func- 
tion, £(r), on scales 10 to 100 Mpc, having £(r) too little power on the large 
scales compared to the power on smaller scales. The APM survey (Maddox 
et al. 1990, 1991), giving the galaxy angular correlation function, the 1.2 
Jy IRAS power spectrum, the QDOT survey (Saunders et al. 1991), X- 
ray observations (Lahav et al. 1989) and radio observations (Peacock 1991; 
Peacock &: Nicholson 1991) agree with the quoted conclusion. As shown in 
recent studies of galaxy clustering on large scales (Maddox et al. 1990; Ef- 
stathiou et al. 1990b; Saunders et al. 1991) the measured rms fluctuations 
within spheres of radius 20/i -1 Mpc have value 2-3 times larger than that 
predicted by the CDM model. 

Alternative models with more large-scale power than CDM have been intro- 
duced in order to solve the latter problem (Peebles 1984; Shah & Stecker 
1984; Valdarnini k Bonometto 1985; Bond et al. 1988; Holtzman 1989; 
Efstathiou et al. 1990a; Schaefer 1991; Turner 1991; Holtzman & Primack 
1993; Shaefer & Shan 1993; White et al. 1993). 

In this work, we show how the discrepancy between the predicted and ob- 
served X-ray temperature distribution function of clusters may be reduced 
when non-radial motions, that develop during the collapse process, are taken 
into account. We compare the distribution function obtained taking into 
account non-radial motions with X-ray temperature distributions obteined 
using a simple CDM spectrum. 

The plan of this paper is the following: in § 2 we find the effects of non- 
radial motions on the X-ray temperature distribution function and then we 
compare this prevision to that by a simple CDM model. In § 3 is devoted 
to the discussion of our results. 



2 X-ray temperature distribution function 

The Press-Schechter theory provides an analytical description of the evolu- 
tion of structures in a hierarchical Universe. In this model the linear density 
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field, p(x,t), is an isotropic random Gaussian field, the non-linear clumps 
are identified as over-densities (having a density contrast 5 C ~ 1.68; Gunn & 
Gott 1972) in the linear density field, while a mass element is incorporated 
into a non- linear object of mass M when the density field smoothed with a 
top-hat filter of radius Rf, exceeds a threshold 8 C (M oc R^). The fraction 
of mass in objects of mass greater than M is given by: 

while the comoving number density of non-linear objects of mass M to 
M + dM is given simply by the derivative of F(> M): 

where p is the mean mass density, v = -^jjjj and a(M) is the linear power 
spectrum evaluated at the epoch when the mass function is desidered, given 
by 

1 f°° 

a 2 (M) = — - / dkk 2 P{k)W 2 (kR) (3) 
2tt z Jo 

where P(k) is the power spectrum and W(kR) is a top-hat smoothing func- 
tion: 

W(kR) = -j^3 (sin kR - kR cos kR) (4) 

The redshift dependence of Eq. (2) can be obtained remembering that 

Sc(z)D(0) 



a (M)D(z) 



(5) 



being D(z) the growth factor of the density perturbation. In Eq. (2) 
Press-Schechter introduced arbitrarily a factor of two because 

J™dF{M)= l - (6) 

so that only half of the mass in the Universe is accounted for. Several 
papers have tried to solve the problems of this model in order to improve it 
(Colafrancesco et al. 1989; Peacock & Heavens 1990; Bower 1991). In spite 
of the quoted problems, N-body simulations developed by Efstathiou et al. 
(1988), Brainerd & Villumsen (1992) showed that Press-Schechter analytic 
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theory provides a good description of the evolution of the distribution of 
mass amongst groups and clusters of galaxies (multiplicity function) and this 
has proven particularly useful in analyzing the number counts and redshift 
distributions for QSOs (Efstathiou & Rees 1988), Lyman a clouds (Bond et 
al. 1988) and X-ray clusters (Cavaliere & Colafrancesco 1988). 
The mass function can be transformed into the temperature distribution if a 
relation between the mass and the temperature is known. For a flat matter- 
dominated Universe a simple scaling gives (Kaiser 1986; Evrard 1990): 

io^) {1 + z) (7) 

The mass variance present in Eq. (2) can be obtained remembering the 
Eq. (3) and once a power spectrum, P(k) = Ak n T 2 (k) is fixed giving the 
transfer function T(k) : 

T(k) = [In (1 + 4.164A;)] 2 • (192.9 + 1340/c + 

+ 1.599 • 10 5 /c 2 + 1.78 • W 5 k 3 + 3.995 • 10 6 fc 4 ) _1/2 (8) 

(Ryden & Gunn 1987) and A is the normalizing constant. The normalization 
is obtained, as usual, imposing that the variance in a sphere of 8h~ 1 Mpc, erg, 
is unity. As shown by Bartlett & Silk (1993) the X-ray distribution function 
obtained using a standard CDM spectrum over-produces the clusters abun- 
dances. The discrepancy can be reduced taking into account the non-radial 
motions that originate when a cluster reaches the non-linear regime. In fact, 
the Press-Schechter temperature distribution requires specification of 5 C and 
the temperature-mass relation. The presence of non-radial motions changes 
both S c and the T-M relation. As shown by Barrow & Silk (1981) and Szalay 
& Silk (1983) the gravitational interaction of the irregular mass distribution 
of proto-cluster with the neighbouring proto-structures gives rise to non- 
radial motions, within the protocluster, which are expected to slow the rate 
of growth of the density contrast and to delay or suppress the collapse. 
According to Davis & Peebles (1977) the kinetic energy of the resulting 
non-radial motions at the epoch of maximum expansion increases so much 
to oppose the recollapse of the proto-structure. Numerical N-body simula- 
tions by Villumsen & Davis (1986) showed a tendency to reproduce this so 
called previrialization effect. In a more recent paper by Peebles (1990) the 
slowdown of the growth of density fluctuations and the collapse suppression 
after the epoch of the maximum expansion were re-obtained using a numer- 
ical action method. In the central regions of a density peak (r < O.bRf) 
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the velocity dispersion attain nearly the same value (Antonuccio-Delogu & 
Colafrancesco 1995) while at larger radii (r > Rf) the radial component is 
lower than the tangential component. This means that motions in the outer 
regions are predominantly non-radial and in these regions the fate of the 
infalling material could be influenced by the amount of tangential velocity 
relative to the radial one. This can be shown writing the equation of motion 
of a spherically symmetric mass distribution with density n(r): 

±L n{Vr) + d_ n ( v *) + (2{ v l) - ( V 2)) Vi + n(r)^(v r ) = (9) 

where (ty) and (v$) are, respectively, the mean radial and tangential stream- 
ing velocity. Eq. (9) shows that high tangential velocity dispersion ({v^} > 
2(v 2 )) may alter the infall pattern. The expected delay in the collapse of a 
perturbation may be calculated solving the equation for the radial acceler- 
ation (Peebles 1993): 

if-w-^ (10) 

where L(r, v) is the angular momentum and g(r) the acceleration. Writing 
the proper radius of a shell in terms of the expansion parameter, a(ri,t): 

r(n,t) =ria(n,t) (11) 

remembering that M = ^-pf,(ri,t)a 3 (ri,t)rf, that ||^ = pb, where H is the 
Hubble constant, and assuming that no shell crossing occurs so that the 
total mass inside each shell remains constant, (p(ri,t) = a^'t) ) Eq- (10) 
may be written as: 

fa H 2 {l + 5) AG 2 L 2 

dt 2 ~ 2a 2 + H*(l +3) 2 r} a 3 ( ' 

or integrating the equation once more: 



( da 2 _ 2 



~l + 5~ 




+/ 


a 





H 



_ 2 -^da-2C (13) 



M° (1 + 7)' 



a" 



where 5 = Pt po po is the overdensity and C is the binding energy of the shell. 
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Figure 1: The threshold <5 C in function of the mass M, for a CDM spectrum 
with Rf = 3Mpc, taking account of non-radial motions. 
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Using a technique similar to that by Bartlett & Silk (1993) it is possible 
to obtain the overdensity at the time of collapse: 



1 + 



8G 2 



figif 6f4°$(l + <S) 2 



Qmax £ 2 . ^Q, 







a;' 



(14) 



where <5 co = 1.68 is the critical threshold for a spherical model. Filtering 
the spectrum on cluster scales, Rf = ?>h~ l Mpc, we obtained the total spe- 
cific angular momentum, h(r,u) = L(r,v)/M s h, acquired during expansion, 
integrating the torque over time (Ryden 1988a, Eq. 35): 



h(r, v) 



r t b 



5/2 



(1-COS0)' 



48M sh Jo (0 - sin tf) 4/3 h W - h (0) | 



(15) 



The functions /i(i?), /2(^) ar e given by Ryden (1988a) (Eq. 31), while the 
mean overdensity inside the shell, 5(r), is given by Ryden (1988b): 



S(r,v) = -r / duo - (5(cr) 



(16) 



The mass dependence of the threshold parameter, 8 c (y), can be found as 
follows: we calculate the binding radius, r&, of the shell using Hoffmann & 
Shaham's criterion (1985): 

T c {r, v) < t (17) 

where T c (r, v) is the calculated time of collapse of a shell and t Q is the Hubble 
time. We find a relation between v and M through the equation M = 
■iirprfj?). We so obtain 5 c (v(M)). In Fig. 1 we show the variation of the 
threshold parameter, 5 C (M), with the mass M. Non-radial motions influence 
the value of 5 C increasing its value for peaks of low mass while leaving its 
value unchanged for high mass peaks. As a consequence, the structure 
formation by low mass peaks is inhibited. In other words, in agreement with 
the cooperative galaxy formation theory (Bower 1993), structures form more 
easily in over-populated regions. To get the temperature distribution it is 
necessary to know the temperature-mass relation. This can be obtained 
using the virial theorem and energy conservation (Bartlett &; Silk 1993). 
Adopting the same notation of Bartlett & Silk we have that: 



T = (6AkeV) 



M ■ h 

1O 15 M 



2/3 



1 + 



v^J* L2dr 



M 2 r* 



(G 2 ff 2 fi M 2 )i/3 



(18) 
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where rj = r ta /xi, being r ta the radius of the turn-around, while x\ is de- 
fined by the relation M = 47rxf / o/3 and ip = r e ff/r ta where r e ff is the 
time-averaged radius of a mass shell. Eq. (18) was normalised to agree with 
Evrard's (1990) simulations for L = 0. 

In Fig. 2 the temperature distribution derived taking into account non-radial 
motions is compared with the X-ray temperature distribution predicted by 
a simple CDM model. As shown the remarkable discrepancy between pure 




Figure 2: X-ray temperature distribution function. The dashed line gives the 
temperature function for a pure CDM model, with Rf = 3Mpc. The solid 
line is the same distribution but now taking account of non-radial motions. 

CDM previsions and a CDM that taking account the non-radial motions is 
evident. In a our paper (Del Popolo h Gambera 1997) we compare our X- 
ray temperature distribution function with the experimental data obteined 
by Edge et al. (1990) and by Henry & Arnaud (1991). 
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3 Conclusion 



In this paper we have shown how non-radial motions may reduce the dis- 
crepancy between the observed X-ray temperature distribution function of 
clusters and that predicted by CDM model. As shown even with low normal- 
ization, CDM power spectrum over-produces the cluster abundances. The 
non-radial motions present in the outskirts of protoclusters have the effect to 
reduce the cluster abundances producing a reasonable cluster temperature 
distribution. Another feature of this model that shall be studied in a next 
paper is the role of non-radial motions on the two-point correlation function 
on scales larger than WMpc (Gambera & Del Popolo in preparation). We 
expect that our model shall be able to predict acceptable correlation func- 
tions. 
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